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^§ : Abstract 

Based on the Renormalization Group method, a reduction of non integrable 
multi-dimensional hamiltonian systems has been performed. The evolution equa- 
tions for the slowly varying part of the angle-averaged phase space density, and 
for the amplitudes of the angular modes have been derived. It has been shown 
that these equations are precisely the Renormalization Group equations. As an 
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application of the approach developed, the modulational diffusion in one-and-a 
half degree of freedom dynamical system has been studied in detail. 
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1. Introduction 

It is well-known that dynamical systems may exhibit irregular motion in certain regions of 
phase spaceEra. These regions differ in size, from being considerably small, to occupying 
large parts of phase space. This depends mostly on the strength of the perturbation, as well 
as on the intrinsic characteristics of the system. For comparatively small perturbations the 
regularity of the motion is expressed in the existence of adiabatic action invariants. In the 
course of nonlinear interaction the action invariants vary within a certain range, prescribed 
by the integrals of motion (if such exist). For chaotic systems some (or all) of the integrals of 
motion are destroyed, causing specific trajectories to become extremely complicated. These 
trajectories look random in their behavior, therefore it is natural to explore the statistical 
properties of chaotic dynamical systems. 

Much experimental and theoretical evidenceEm of nonlinear effects in the dynamics of 
particles in accelerators and storage rings is available at present. An individual particle 
propagating in an accelerator experiences growth of amplitude of betatron oscillations in a 
plane transverse to the particle orbit, whenever a perturbing force acts on it. This force may 
be of various origin, for instance high order multipole magnetic field errors, space charge 
forces, beam-beam interaction force, power supply ripple or other external and collective 
forces. Therefore, the Hamiltonian governing the motion of a beam particle is far from 
being integrable, and an irregular behavior of the beam is clearly observed, especially for a 
large number of turns. 

The idea to treat the evolution of chaotic dynamical systems in a statistical sense is 
not new; many rigorous results related to the statistical properties of such systems can be 
found in the book by Arnold and Aveza. Many of the details concerning the transport 
phenomena in the space of adiabatic action invariants only are also well understoodn. In 
this aspect the results presented here are in a sense re-derivation of previously obtained ones 
by means of a different method. What is new however, is the approach followed to obtain 
the diffusion properties in action variable space, as well as a new evolution equation for the 
angle-dependent part of the phase space density. Furthermore, instead of the widely used 
phenomenological method to derive the diffusion coefficient (tensor), the procedure pursued 
in the present paper is more consistent one, with a starting point the Liouville's equation 
for the phase space density. 

We first employ the Projection Operator method of Zwanzigu to derive the equations for 
the two parts of the phase space density: the averaged over the angle variables part F, and 
the remainder G [see Eq. ( j2.10| ) in the next Section]. As expected, the two equations are 
coupled. Next we extract the relevant long-time scale behavior embedded in the equations 
for F and G by means of the Renormalization Group (RG) methodQ'0. It is remarkable, and 
at the same time not surprising that the equations governing the long-time scale dynamics 
are the Renormalization Group equations (RGEs). These are obtained in Section 4 through 
renormalizing the perturbative solution of the equations for F and G [see Eqs. fl2.22|) and 



( |2.2^ ) of Section 2]. Finally, in Section 5 one-dimensional example of a chaotic system is 
considered to demonstrate the approach developed here. 



2. Projection Operator Method 

Single particle dynamics in cyclic accelerators and storage rings is most properly described by 
the adiabatic action invariants (Cour ant- Snyder invariantsu) and the canonically conjugate 
to them angle variables. However, to be more general we consider here a dynamical system 
with N degrees of freedom, governed by the Hamiltonian written in action-angle variables 

(J, ex.) as 

H(ex,J;9)=H (J) + eV(ex,J;9), (2.1) 

where 9 is the independent azimuthal variable (widely used in accelerator physics), playing 
the role of time and J and ex are iV-dimensional vectors 

J= (Ji,J2,---,Jn), a = (ai,a 2 ,...,a N ). (2.2) 

Moreover H (J) is the integrable part of the Hamiltonian, e is a formal small parameter, 
while V(a, J; 9) is the perturbation periodic in the angle variables 

V{a, J; 9) = £V m (J; 0) exp (im ■ a), (2.3) 

m 

where J2' denotes exclusion of the harmonic m = (0, 0, . . . , 0) from the above sum. The 
Hamilton's equations of motion are 

da k dV dJ k dV 

- W = -o k (J) + e— , ^T = - e ^? (2 - 4) 

where 

U 0k (J) = jr^. (2.5) 

dJk 

In what follows (in particular in Section 4) we assume that the nonlinearity coefficients 

^ (j > - St, < 2 - 6) 

are small and can be neglected. The Liouville's equation governing the evolution of the 
phase space density P(ct, J; 9) can be written as 

^-P(a,J;0)= \£o + e£ v (e)]P{a,J;e). (2.7) 

o9 L J 

Here the operators Cq and C v are given by the expressions 

r d ? dV d dV d 

da k oa k dJ k dJ k da k 

where summation over repeated indices is implied. Next we define the projection operator 
onto the subspace of action variables by the following integral: 



2tt 2tt 



^( J ! °) = TT^n / dai . . . / dajv/(a, J; 0), (2.9) 

[Z7T) o o 

where /(a, «/; 0) is a generic function of its arguments. Let us introduce also the functions 
F = VP, G=(l-fi)P = CP, (P = F + G). (2.10) 

From Eq. (|2.7|) with the obvious relations 

VjCo = jC V = 0, VjC v P = (2.11) 

in hand it is straightforward to obtain the equations 

f= £ ^ G =^E4 (212) 

BC 

— = C G + eCC v G + eC v F (2.13) 



Our goal in the subsequent exposition is to analyze Eqs. Q2.12Q and ( |2.13|) using the 



RG method. It will prove efficient to eliminate the dependence on the angle variables in G 
and V by noting that the eigenfunctions of the operator jCq form a complete set, so that 
every function periodic in the angle variables can be expanded in this basis. Using Dirac's 
"bra-ket" notation we write 

\ n ) = / 27r W/2 eX P( in " a )' H = / 2 xjV/2 eX P(~ in " a )- ( 2A4 ) 

The projection operator V can be represented in the formBj 

p = p = |0)(0|. (2.15) 

One can also define a set of projection operators V n according to the expression!^ 

V n =\n)(n\. (2.16) 

It is easy to check the completeness relation 

E^« = l, (2.17) 



from which and from Eq. ( |2.15p it follows that 



C= $3'|n)(n|. (2.18) 



Decomposing the quantities F, G and V in the basis ( ^.14[ ) as 



F = F(J;6)\0), (2.19) 



G(a,J;9)= £ 'G m (J;fl)|m), 



(2.20) 



V(a,J;9) = J^'v n (J;9)\n), 

n=£0 



(2.21) 



from Eqs. ( |2.12| ) and ( |2.13| ) we obtain 

dF 
86 



dF . d /w 



<9J fc 



(2.22) 



06 



-in k uj ok G n + ie^ 



n k V n 



dGr, 

1 dJu 



d 
m k Trr(yi 



dJ k 



n—m K -' m 



G r 



+ ien k V n 



dF 

dJ k 



(2.23) 



The above equations comprise the the starting point in the analysis outlined in Section 
4. We are primarily interested in the long-time evolution of the original system governed by 
certain amplitude equations. These will turn out to be precisely the RG equations. 



3. Renormalization Group Reduction of Hamilton's Equations 



Let us consider the solution of Hamilton's equations of motion (|2.4j ) for small perturbation 
parameter e. It is natural to introduce the perturbation expansion 



(0) , (1) , 2 (2) . 

(Xk — oi k + ea k + e a k + 



Jk - J k + eJ k + e J k + 



The lowest order perturbation equations have the trivial solution: 



a k = u ok 9 + Lp k , 



T {0) - Au 



(3.1) 



(3.2) 



where A k and (p k are constant amplitude and phase, respectively. We write the first order 
perturbation equations as 



da 



(i) 



dO 



im(A)4 1) + 



dV 
dT k 



dJ, 



(i) 



dV 



d9 " da f 



(3.3) 



Assuming that the modes V n ( J; 9) are periodic in 9, we can expand them in a Fourier series 

oo 

= J2 V n (J;fi)exp(itJtv n 0). (3.4) 



V n (J; 



/i= — oo 



If the original system ( |2.1| ) is far from primary resonances of the form: 
n k 'u ok + fxu R = 



(3.5) 



we can solve the first order perturbation equations ( |3~3"D , yielding the result: 



a k = % 1^ 2JluyA)miV m {[i) ■ -2 exp (im a (p a ) 

m 11 {m s u 0s + fJ,U m ) 

- <E / E g # j ^ [l{msU ° S + ^ m)9] ew <im.<p.), (3.6) 

^y dA k m 8 Uo 8 + tii/ m 

4" = -tt>H,VM 6XP [ ' (m ' a, °: + "^ exp (im, V ,), (3.7) 

The second order perturbation equations have the form: 



' 2) ^ (A\ T^ + X ^ 7fc/ f« rW + ^ rW + g2y n .(i) fQ s^ 



d£ "" v 7 ' 29i 8 ' s dA fc dA. ' M fe aa{ 0) 

d4! _ ^ (1) _ gy (1) 



The solution to Eq. ( |3.9|) reads as 

+ 2vrft(#) £ '£m fe m*m s7 * s (A)|K„Gu)| 2 — 3fc( 7 ; a) 
+ oscillating terms, (3.10) 

where 

vr3?(x; j/) = ~^— 2 , limK(x; j/) = <J(y), (3.12) 



and 7 is a small real quantity added ad hoc in the denominators of the expressions (|3.6j ) and 
( p.7|) . The limit 7 — > will be taken in the final result. 

As expected, in the second order perturbation solution (|3.10p the first and the second 



terms are secular, because 1Z(6)= 9. To remove these secularities we follow the general 
prescription of the RG methodQS. First, we select the slowly varying part of the perturba- 
tion solution governing the long-time evolution of the system. Up to second order in the 
perturbation parameter e it consists of the constant zero order term A^ and the second 
order secular terms. Next, we introduce the intermediate time r, and in order to absorb the 
difference r = 9 — (9 — r) we make the new renormalized amplitude A^r) dependent on 
r. Since the long-time solution thus constructed should not depend on r its derivative with 
respect to r must be equal to zero. This also holds for r = 9, so that finally 



dA k (6) 2 , d\V m (A; fi)\ 2 „, , 

— — — = 2vre > 2Jm k mi — 3fc(7; m s w 0s + /x^ m ) 

Clt7 „.^ n ,, OAl 



m>0 M 
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27rg2 H Y^ m k m i m slis (A) \V m (A; n)] 2 ^-^(7; 



(3.13) 



a=m r UOr+A«'m 



Equation (|3.13| ) is known as the Renormalization Group equation (RGE). It describes 
the slow long-time evolution of the action variables. 



4. Renormalization Group Reduction of Liouville's Equation 



We consider the solution of Eqs. (|2.22|) and (|2.23j ) for small e by means of the RG method. 
For that purpose we perform again the naive perturbation expansion 



F = F<°> + eF« + e 2 F^ + ■■■ 



G n = G® + eG« + e 2 ^ + ■■■ 



(4.1) 



and substitute it in Eqs. Q2.22|) and (|2.23|) . The lowest order perturbation equations have 
the obvious solution 



F<® = F (J), G^ = W n {J) exp (-in k oo ok t 

The first order perturbation equations read as: 

J2 n k V n W-nexp (iniUoiS 



dF™ d 

1 



oe 



dJ k 



(4.2) 



(4.3) 



dGU 

dd 



-in k u ok 



Gi 1 ' + in k V T 



8J k 



+C 



n k V n 



dWr, 

1 dJ k 



d (V 



dJ k 



n—m'Vm) 



exp {-imiujoit 



(4.4) 



We again assume that the modes V n (J; 9) are periodic in 9, so that they can be expanded 
in a Fourier series ( |3.4| ). If the original system ( |2.1D exhibits primary resonances of the form 
( p.5|) in the case when uqu does not depend on the action variables, we can solve the first 
order perturbation equations (|4.3|) and fl4.4|) . The result is as follows: 



F (1) = iK{6) J2 



" (R) 
n 



i(«) 



d_ 
dJk 



V n (R) 



(R) 



W_ n (R) 



d 



+ T!T!nk Wk [vMW- 



n fj, 



exp[i(nitJoi + fiu n )9] 
riiuoi + \xv n 



(4.5) 



G$ = £ n exp (-in k uj 0k t 
where 



(4.6) 



r "d(q\x dF °if ( n "a 
G n = tK.{9)d nn (R)n k ——V n \ 



+ m{B) 



E 

n — n (R) 



" ' TA I n i 



(R) 



UQl \ 



dW„ 



»(«) 



vr j 



dJ k 



f (R)\ d 



(R) 

V„iO) | | W n _ n n ! 

VR 



9fow i; , v exp [i(n/o;o« + /^n)0] 



d.J k 



niUJoi + yuz/ r 



+ LL "fcK*-™^)^ m k —-[V n - m (fj)W m ] \ r (4.7) 

1 oJ k dJ k I (m - mijujoi + ^v n -m 



m [i 



In the above expressions J2" denotes summation over all primary resonances ( |3.5| ). To obtain 
the desired RG equations we proceed in the same way as in the previous Section. The first 
order RG equations are 



36 



E 



ze > n 



dJ k 



[V n(R) {J; n R )W_ n(R ){J)}, 



(4.8) 



-gg- = K°nn&) n k-Qj-V n {J; n R ) 



+ ze £ "{^^(J;/^) ^-^' - (n fc - n k R) )^-[V n{R) {J^R)W n _ n(R) {J)]\, (4.9) 



where 



Mb 



(R) 

vr 



(4.10) 



Equations ( |4.8| ) and ( |4.9| ) describe the resonant mode coupling when strong primary reso- 
nances are present in the original system. 

Let us now assume that the original system is far from resonances. Solving the second 
order perturbation equation for F^ and G^} 



dF^ . d ( , (1) 



(4.11) 



dG<® 



dF (D 



00 
we obtain 



+ iUkUokG® = in k V n —— + iJ2 



dJ k 



n k V 7 



dG£) 







dJk dJi V 



V G (1) 



(4.12) 



F (2) = 2vr7e(#) 



g 

dJ k 






n>0 A 



+ oscillating terms, (4.13) 



^ = ^T»exp (-m s u; 5 



(4.14) 



where 

V (u) <9 2 

^ = in{B)n k n{£ "^ \y*Qi)W n 



+ *W)E'E 



•£ [n s - m s )u 0s + [iv n -m \ oJk \ oJi 



pp. F) ( f)W 

+ n k niV n - m (ix) (Cm(i")^r.) + m k m l —l\V n _ m (fi)\ 7 



i2 uyv n 



dJi 



-m k ni 



d 
dJ k 



d_ 
'dJi 



v n -M—(v:_ m (fi)w n 



+ oscillating terms, 



(4.15) 



and the functions 11(9) and 3ft(x; y) are given by Eqs. (|3.11|) and (|3.12|) , respectively. It is 
now straightforward to write the second order Renormalization Group equations. They are: 



89 dJ k 



E E n ^lKi(A)| 2 ^9?(7; n s u 0s + Xi/ n ) 
k Ln>0 a aj i 



(4.16) 



dW, 



39 



n ■ 2 

- ie 



V™/E 



n s u 0s + fiu n 8J k dJi 



\v:^)w n 



+ - 2 E'E- 



1 



m fi 



n s - m s )uj 0s + [iv r 



-n k miV n ^ m ([i)—\V*_ m (fi)-^j- 



pp. p. I FiXAl 

+ n k riiV n _ m (fi) - - - - (y*_ m (fj,)Wr^ + m k mi——i \V n _ m (^)\ 2 



,2 uv ' n 



dJkdJl 



W-Jl 



dJi 



-m k ni 



d 
dJ k 



Vr^(ji)^(VLMW T 



(4.17) 



The RG equation ( |4. 16|) is a Fokker-Planck equation describing the diffusion of the 
adiabatic action invariant. It has been derived previously by many authors (see e.g. the 
book by Lichtenberg and Liebermana and the references therein). It is important to note that 
our derivation does not require the initial assumption concerning the fast stochastization 
of the angle variable. The fact that the latter is indeed a stochastic variable is clearly 
visible from the second RG equation (|4.17|) , governing the slow amplitude evolution of the 
angle-dependent part of the phase space density. Nevertheless it looks complicated, its most 
important feature is that equations for the amplitudes of different modes are decoupled. In 
the case of isolated nonlinear resonance Eq. ([4.17[) acquires a very simple form as will be 



shown in the next Section. 



5. Modulational Diffusion 



As an example to demonstrate the theory developed in previous sections, we consider the 
simplest example of one-and-a-half degree of freedom dynamical system exhibiting chaotic 
motion 

H (J) = XJ + H 8 (J), V (a, J; 9) = V (J) cos (a + £ sin vd). (5.1) 



The Hamiltonian (|5.1| ), written in resonant canonical variables describes an isolated nonlin- 
ear resonance of one-dimensional betatron motion of particles in an accelerator with mod- 
ulated resonant phase (or modulated linear betatron tune). The modulation may come 
from various sources: ripple in the power supply of quadrupole magnets, synchro-betatron 
coupling or ground motion. The resonance detuning A defines the distance from the reso- 
nance, S is the amplitude of modulation of the linear betatron tune and v is the modulation 
frequency, where £ = H/V. Without loss of generality we consider £ positive. Since 

^ = A + ^T' V 1 (Jm) = ±V(J)J lt {t), (5.2) 



where J n {. z ) is the Bessel function of order n, the RG equation ( |3 . 1 3| ) for the amplitude A 

can be rewritten as 



,u '"" ' ° W'wjud ' ' ' '' 



d6 2v dA 



1{A)V\A)-[J>{0 



(5.3) 



Here the square brackets [z] encountered in the index of the Bessel function imply integer 
part of z. Moreover, in deriving the expression for Vi( J; fi) in Eq. ( |5.2| ) use has been made 
of the identity 

oo 

exp (iq sin z) == Yl X(|<?|) exp [mzsgn(g)], (5.4) 

n=— oo 

and finally, the limit 7 — > in Eq. (|3.13|) has been taken. For small value of £ utilizing the 
approximate expression for the derivative of Bessel functions with respect to the order we 
obtain 
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cL4 

~d6 



2v \dA 



V\A)J?^(0 



2J^ ] (0\n(^j 1 (A)V\A) 



(5.5) 



Let us now turn to the RG equations ( |4.16| ) and ( |4.17| ). They can be rewritten in the 
form: 



8F Tie 2 d 
~~cW~ ~2^7dJ 



v\J)J^(0 



dJ 



(5.6) 



1 dW n 
W„ 06 



me 2 n 



, ,^d_/ dV^ , ne 2 n ^ 2 , d f^dV 



2v sin (ttlu / vj J -^^d J \ dJ ) + 2v 



■ J m^ij[ v: A j 



(5.7) 



Equation (|5.7|) suggests that the amplitudes W n of the angular modes G n exhibit exponential 
growth with an increment 



r -£^]«>£^ 



(5.8) 



Equation (|5.6|) is a Fokker-Planck equation for the angle-independent part of the phase space 
density with a diffusion coefficient 



V{J) = 7 ^;V 2 {J)J 2 m {fi). 

In Figures [I|-|3] the reduced diffusion coefficient 

22 



V^ J, 



v) ?re 2 V 2 (J) 



V(J) 



(5.9) 



(5.10) 



has been plotted as a function of the ratio between the amplitude and the frequency of 
the modulation. Three typical regimes corresponding to different values of A/S used as a 
control parameter have been chosen. In the first one depicted in Figure |l] we have taken 
the resonance detuning twice as large as the amplitude of the modulation (A = 22). In this 
case there is no crossing of the main resonance described by the Hamiltonian ( |5.1|) and the 
diffusion coefficient decreases very rapidly after reaching its maximum value at £ = 0.25. 
The cases of periodic resonance crossings for A = 5 and A = 2/2 are shown in Figure ^| and 
Figure |3], respectively. 



6. Concluding Remarks 

In the present paper we apply the Renormalization Group method to the reduction of non 
integrable multi-dimensional Hamiltonian systems. The notion of reduction is used here in 
the sense of slow, long-time behavior, that survives appropriate averaging and/or factorizing 
of rapidly oscillating contributions to the dynamics of the system. 

As a result of the investigation performed we have derived evolution equations for the 
slowly varying part of the angle-averaged phase space density, and for the amplitudes of 



11 



the angular modes. It has been shown that these equations are the Renormalization Group 
equations. 

The case of a one-dimensional isolated nonlinear resonance with a resonant phase (or 
linear unperturbed tune) subjected to periodic modulation has been studied in detail. The 
coefficient of modulational diffusion, as well as the exponential growth increment of the 
amplitudes of angular modes have been obtained in explicit form. 
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Fig. 1. Reduced diffusion coefficient T>^ R > as a function of the ratio £ = H/V for A = 2; 
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Fig. 2. Reduced diffusion coefficient T>^ > as a function of the ratio £ = z>/v for A = H. 
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Fig. 3. Reduced diffusion coefficient T>^ > as a function of the ratio £ = 3/v for A = S/2. 
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